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Abstract. We have used the Westerbork Synthesis Radio Telescope to detect HI 21cm absorption at z ~ 0.7645 in the grav- 
itational lens system towards PMN J0134-0931. The 21cm profile has two broad components, with peak optical depths of 
0.047 ± 0.007 and 0.039 ± 0.007, at heliocentric redshifts 0.76470 ± 0.00006 and 0.76348 ± 0.00006, respectively. The redshift 
of the stronger component matches that of Call H and K absorption detected earlier. The absorption has a total velocity width 
of ~ 500 km s"' (between nulls) and an equivalent width of 7.1 + 0.08 km s"'. This would imply a total HI column density of 
2.6 + 0.3 X 10^' cm"-, for a spin temperature of 200 K and a covering factor of unity. The high estimated HI column density 
is consistent with the presence of large amounts of dust at the lens redshift; the intervening dust could be responsible for the 
extremely red colour of the background quasar. 
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1. Introduction 

The z ~ 2.2 radio-loud quasar PMN J0134-0931 is an ex- 
ceedingly interesting object. Its radio continuum consists of 
at least 6 components, with a maximum image separation of 
~ 700 mas f Winn et al. 20Q2'(W02):' Winn et al. 2003I (W03)). 
Five of these images (A - E) have the same spectral index 
and a curved arc of emission is seen between A and B. These 
properties caused W02 to argue that the object is a gravi- 
tationally lensed system. Similarly, Gregg et al. (2002) used 
the multiplicity of components in a near-IR image as well as 
the high inferred source luminosity as evidence for lensing. 
Interestingly enough, PMN JO 1 34-093 1 is only the second lens 
of galactic mass to show an image multiplicity greater than 
four; the other such lens is the six-component system towards 
CLASS B1359-i-154(Rusin et al. 2001). 

Call H and K Unes have been detected towards 
PMN JO 1 34-093 1 , at z = 0.7645 1 + 0.000 1 6 (' Hall et al. 1 (Mk 
the lens is thus likely to be at this redshift. The Call absorption 
is very strong, with rest frame equivalent width ~ 5.3 A, sim- 
ilar to that seen in the lens B0218-I-357 (Browne et al. 1993). 
Further, the quasar is highly reddened in optical wavebands, 
with B-K > 1 1 ( Greg g et al. 2002^ ; this suggests the presence 
of large amounts of dust, either in the host galaxy of the quasar 
or in the lensing galaxy. While the differential extinction of the 
quasar images suggests that at least some of the reddening is 
due to dust in the absorber, it has so far not been possible to 
conclusively determine the location of the dust. Assuming the 
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dust is at z = 0.7645 implies E(B - V) - 1.315, i.e. a visual 
extinctionAv = 3.85 (Hall et al. 2002.1. In conjunction with the 
strong Call lines, this would suggest that the lens is a gas-rich, 
late-type galaxy. 

More recently, Keeton & Winn (2003) have come up with 
the first quantitative model of the lensing system that explains 
most of its observed properties. The background source is as- 
sumed to have two components and the lensing system to con- 
sist of two galaxies. The observed 5 components A - E are 
then the images of one of the background sources and com- 
ponent F is an image of the second, with the remaining two 
images of the latter source too faint for detection. The model 
is supported by new high-resolution radio and optical observa- 
tions (W03), which found that F has a different spectral index 
from images A - E. These observations also showed that the 
observed large sizes of components C and E are not intrinsic but 
instead arise due to scatter broadening at the lower frequencies. 
Finally, W03 find marginal evidence for a detection of the two 
lens galaxies in Hubble Space Telescope (HST) images of the 
system; Keeton & Winn (2003) predict these to be spiral galax- 
ies, each with a velocity dispersion of about ~ 120 km s"', and 
a projected separation of ~ 0.4". 

We present, in this Letter, Westerbork Synthesis Radio 
Telescope (WSRT) observations of HI 21cm absorption at z ~ 
0.7645 towards PMN J0134-0931. The usual association of 
dust with high column densities of neutral gas confirms that 
the dust is likely to be concentrated in the intervening lens sys- 
tem. We describe the WSRT observations and data analysis in 
section|3 the spectra, results and models are presented and dis- 
cussed in section|5] 
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2. Observations and data analysis 

PMN JO 134-0931 was observed with the UHF-High band of 
the WSRT on March 28 and 31 and July 5 and 7, 2003, with 
the DZB correlator as the backend. A single 5 MHz IF band 
was used on all occasions, sub-divided into 1024 channels and 
centred at an observing frequency of 805 MHz; this yielded a 
frequency resolution of ~ 4.88 kHz, i.e. a velocity resolution of 
~ 1.8 km s"'. The standard calibrators 3C48 and 3C147 were 
used to calibrate the absolute flux density scale at all observing 
epochs; their adopted flux densities are on the scale of Baars et 
al. (1977). The system bandpass was calibrated with observa- 
tions of 3C48, 3C147, 3C454.3 or 3C138. All fourteen WSRT 
antennas were available for the observations in all observing 
runs. The total on-source time was 17.6 hours. 

The data were converted from the telescope format to FITS 
and then analyzed in classic AIPS, using standard techniques. 
Data from the March and July runs were analyzed separately 
and self-calibrated independently. Spectra from the two epochs 
were corrected to the heliocentric frame outside AIPS and then 
averaged together The flux density of PMN JO 134-0931 was 
measured to be 590 + 2 mJy on all observing runs. 

We note that two strong, narrow RFI spikes were found to 
be present not far from the expected location of the redshifted 
21cm line; these were seen in both the March and July runs, 
at the same observing frequencies (i.e. shifted heliocentric fre- 
quencies). The high initial spectral resolution of the observa- 
tions limited the effect of these spikes on neighbouring chan- 
nels; the data were also Hanning smoothed in frequency, to fur- 
flier reduce spectral ringing. 
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Fig. 1. Final WSRT HI 21cm spectrum towards 
PMN JO 134-0931. The x and y axes are heliocentric fre- 
quency (in MHz) and flux density (in mJy) respectively; 
the top axis is heliocentric velocity (in km s '), relative to 
z - 0.7645. The spectrum has been boxcar smoothed by 3 
channels and has a velocity resolution of ~ 10.8 km s ' and an 
RMS noise of ~ 4.1 mJy. 



3. Discussion 

3.1. The 21cm spectrum and derived quantities 

The final WSRT HI 21cm spectrum towards 
PMN JO 134-0931 is shown in figure [T] this has been 
boxcar smoothed by three channels (and re-sampled) and has 
a velocity resolution of ~ 10.8 km s ' and an RMS noise of 
~ 4.1 mJy. The redshift of the Call lines JHall et al. 2002.1 is 
shown below the spectrum, with error bars. 

Two deep features can be immediately identified in the 
spectrum of Fig.^ with their peak optical depths at heliocentric 
frequencies of 804.898 + 0.03 and 805.455 + 0.03 MHz, i.e. he- 
liocentric redshifts of 0.76470+0.00006 and 0.76348+0.00006 
respectively. The redshift of the first of these features agrees 
with that of the Call lines, within the errors. This is also the 
deeper of the two HI components, with a peak optical depth 
of 0.047 + 0.007; the second feature has a peak optical depth 
of 0.039 + 0.007. The HI absorption is extremely wide, with a 
total velocity spread (between nulls) of ~ 500 km s This 
is similar to the width of the Call lines, which have o-y = 
220 + 40 km s'^ JHaU et al. 2002t . 

The equivalent width of the 21cm profile is J T2idV - 
7.06 + 0.08 km s This is is higher than that measured in 
either of the other gravitational lenses in which HI absorp- 
tion has so far been detected, B0218+357 (EW = 2.94 ± 
0.02 km s-i; Kanekar et al. 2003) and PKS 1830-21 (EW = 



5.8 + 0.4 km s '; Chengalur et al. 1999). The large equiva- 
lent width implies a very HI column density A^hi = 2.6 + 
0.3 X 10^1 (rs/200K)(l//) cm-2 (e.g. Kanekar & Chengalur 
2003), where / is the covering factor of the absorber and Ts, 
its spin temperature; this is well into the range observed in 
damped Lyman-a systems (e.g. Wolfe et al. 1986). We note 
that a spin temperature of ~ 200 K is typical of spiral galax- 
ies, both locally dBraun & Walterbos 1992a and at low redshift 
(Kanekar & Ch engalur 2003) . Of course, the above estimate of 
A^Hi gives the sum of HI column densities towards all the differ- 
ent source components giving rise to the 21cm absorption; it is 
not possible to measure the column density towards individual 
source components with the available data. 

The extremely high estimated HI column density implies 
a high dust content in the z ~ 0.7645 lens and explains the 
observed reddening in the optical images ( Gre gg et al. 2002| 
W03). Assuming that the molecular hydrogen column den- 
sity is roughly equal to A^hi, the total hydrogen column den- 
sity A^H = ZA'hj + ^Hi ~ 7.8 X 10^' cm"^ (again assuming 
Ts - 200 K and a covering factor of unity). This implies E(B - 
V) ~ 1 .35, using a Galactic dust to gas ratio; further assuming a 
Galactic extinction law (Ry = 3.1; Binney & Merrifield 19*98) 
results in a visual extinction Ay ~ 4.15. Of course, it is likely 
that the 21cm absorption profile (and hence the estimated col- 
umn density) contains contributions from more than one line 
of sight, with variable optical depth and kinematics across the 
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complex lensing image; it is thus not possible to estimate the 
visual extinction towards individual source components. 

3.2. Modelling the absorption 

The lens towards PMN JO 134-0931 creates a complex radio 
image of six image components (A - F), with a maximum an- 
gular separation of ~ 681 mas (W02). Components A and B 
are seen to connect in a section of lensing arc at faint levels 
in the 1.7 GHz VLBA image (W03). Components C and E 
appear to be strongly scatter broadened at frequencies below 
5 GHz and are almost undetected in VLBA images at 1 .7 GHz 
(W02,W03), as they are over-resolved by even the shortest 
baselines of the VLBA. Note, however, that their flux must be 
present in the integral flux density measured in our spectrum at 
805 MHz; further, they provide additional lines of sight through 
the z ~ 0.7645 absorber. Next, while component F has a sig- 
nificantly steeper spectrum than the other five components, its 
total flux density is only 1 1 ± 1 mJy at 1.7 GHz; it is thus quite 
unlikely to have sufficient flux density at 805 MHz to give rise 
to either of the two main 21cm components. Further, while it 
is possible that the some of the wide, weak absorption arises 
against this component, this would require a fairly high opti- 
cal depth (t 0.5). The 21cm absorption is thus most likely 
to arise against the brighter components A, B, C and D (and 
perhaps the arc of radio emission between A and B), each of 
which have flux densities ^100 mJy at 1.7 GHz. We note that 
component A also shows evidence for scatter broadening at fre- 
quencies below 5 GHz while D is the most strongly reddened 
of the source components; on the other hand, B is the bluest of 
all components and shows no evidence for scatter broadening 
in the VLBA images (W03). This suggests that there is rela- 
tively less dust and gas along the line of sight to component 
B; the two main 21cm features thus appear most likely to arise 
from absorption against components A and D. 

We have explored these possibilities with a simple kine- 
matic model using a differentially rotating disk absorber. The 
modelling does not produce a unique result, but rather succeeds 
in demonstrating that a single, massive galaxy could indeed be 
responsible for the 2 1 cm absorption profile. If a more complex, 
multi-galaxy system were to be invoked, there would be am- 
ple free parameters to fit the integral spectrum reported here. 
Further observational progress will require radio interferome- 
try to separate the absorption along different lines of sight. 

The philosophy of modelling the 21cm profile shape is 
based on the assumption that narrow absorption features occur 
in the spectrum when a compact, high-surface brightness image 
of the background radio source "selects" a small patch from an 
extended, smooth velocity field of a rotating galaxy. Thus, the 
two main narrow 21cm features of Fig. ^are expected to be 
associated with specific source components of the background 
source, once VLBI observations are conducted in the 21 cm 
line. The modelling uses the source structure obtained from the 
1 .7 GHz VLBA maps of W03 for the compact radio compo- 
nents, with a scaling in flux to both match the relative flux den- 
sities of the components in the higher frequency bands, and still 
produce the integral flux density measured at 805 MHz by our 



WSRT observations. The more extended components C and E 
had to be supplemented with diffuse gaussian components to 
account for the flux missing from the VLBA map. 

The two models discussed here (and shown in Fig. |3B]) 
are representative of two broad classes which can reproduce 
the observed 21cm profile. It should be emphasized, however, 
that neither class is able to account for all the observations si- 
multaneously, including the scatter broadening and reddening 
of different source components. Both models are axially sym- 
metric, planar systems with rotation curves that are flat at Vrot 
in the outer regions, with a linear rise from the centre to Re- 
The systemic velocity of the galaxy is offset from 805 MHz 
(z = 0.76448) by Voff, and the galaxy centroid is specified 
by the projected positional offsets AX (west), AY (north) rel- 
ative to image D. The HI optical depth profile for each pixel is 
broadened by a velocity dispersion of 0"^, =7 km s"' before 
computing that pixel's contribution to the integral absorption 
spectrum. 

3.2.1. Model I 

Models belonging to class I are motivated by the need to in- 
flict heavy reddening on image D without reddening B, cou- 
pled with accounting for the scatter broadening of A and C. 
However, this model does not account for the scatter broaden- 
ing of image E. The model parameters are: position angle of 
the receding major axis p = -50° (+p, measured east from 
north), inclination to sky plane / - 85°, Viot - 230 km s"', 
Rc = 1.4 kpc, exponential gas scale height h = 0.7 kpc, 
Voff = 35 km s"', and galaxy position at AX - 490 mas, 
AY = 370 mas. In this highly inclined model, the HI col- 
umn density measured perpendicular to the galaxy plane is only 
3xlO'^(rs/100) cm"^ in the region that covers components A 
and C, but it rises to 3xl02°(rs/100) cm"^ ^j. g j^^^j q j^e 
high disk inclination increases the line of sight column densi- 
ties by an order of magnitude. The 21cm spectrum resulting 
from this model is shown in fig. |3A]. Here, the two narrow 
features arise from absorption against A and D; the extended 
weak absorption wing out to — 1-300 km s"' is achieved by tip- 
ping the intervening disk to a high inclination and thickening 
the HI layer with a scale height of ~ 700 pc (in the spirit of 
Briggs et al. (1985) and Prochaska and Wolfe (1997)). 

3.2.2. Model II 

The second, somewhat orthogonal, class of models (of which 
Model II is an example) has an infinitesimally thin disk, ori- 
ented to cast a shadow on images A, B, C, and E, thereby ac- 
counting for the observed scatter broadening of A, C, and E, 
but implying that all would be vulnerable to reddening in con- 
tradiction with B being the bluest and D being the reddest of 
the images. The model parameters are: p = H-52°, / - 72.5°, 
Vrot - 250 km s"', R^ = 0.8 kpc, (gas scale height h - kpc), 
Voff = -15 km s"', and galaxy position at AX - 535 mas, 
AY - 290 mas. This model has an exponential radial decline 
in HI column density, with scale length 1.5 kpc, from a cen- 
tral peak of 1.5xlO~'(7^s/100) cm"^ measured perpendicular to 
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Fig. 2. Disk galaxy models for the absorption profile. Central Panel [B]: The six principal image components (A-F) of the 
1 .7 GHz VLBA image of W03 (grey-scale), with two sets of ellipses drawn to indicate the orientation and inclination of the 
intervening model disk galaxies. The galaxy centre is marked by a + and a x for Models I and II, respectively; the two filled 
squares indicate the mass centroids of the two mass lens model of Keeton & Winn (2003). Left Panel [A]: 21cm spectrum derived 
from Model I (heavy line) superimposed on the observed 21cm profile (dashed line). Right Panel [C]: 21cm spectrum derived 
from Model II (heavy line) superimposed on the observed 21cm profile (dashed line). In both cases, the zero of the velocity scale 
corresponds to 805 MHz (z - 0.76448, heliocentric) . See text for discussion. 



the disk plane. The spectrum resulting from the above orienta- 
tion is shown in fig. |3C]. Here, the two narrow features arise 
against images A and B, while the broad shoulder at positive 
velocities stems from absorption against image C and the radio 
arc connecting A and B. 

4. Summary 

The 21cm absorption line profile at z ~ 0.7645 towards 
PMN JO 134-0931 has two primary components and a broad 
weak trough. It has a total velocity width of ~ 500 km s"' and 
the largest HI equivalent width so far detected in a gravitational 
lens system. The extremely high estimated HI column density 
A^Hi = 2.6 + 0.3 X 10^' (rs/200K)(l//) cm^^ suggests that the 
reddening of the quasar is due to dust in the z ~ 0.7645 lens. 

None of the galaxy models that assume a single disk-like 
absorber have succeeded in accounting for all the observations, 
although the models can roughly reproduce the observed 21cm 
absorption profile. Specifically, Model I in Sect. l3.2l does not 
account for the observed scatter broadening of component E 
while Model II fails to reproduce the fact that D is the reddest 
and B, the bluest, of all images. Of course, the large 21cm ve- 
locity width of ~ 500 km s"' could also arise if the absorption 
occurs in two separate galaxies. This is in agreement with the 
suggestion of Keeton & Winn (2003) that the lensing system 
consists of two spiral systems. Further, it is consistent with the 
HST image of W03 showing marginal detections of two galax- 
ies, the northern one close to images A, C and E and the south- 
ern one near image D. For example, the northern "galaxy" in 
the HST image might cause the main 21cm feature and most of 
the broad, weak absorption, arising against images A, C, E and 
the radio arc, while the southern "galaxy" causes the secondary 
(blue-shifted) feature (and perhaps some of the broad trough), 
through absorption against image D. The present 21cm data do 



not warrant a detailed modelling of such a two-galaxy absorp- 
tion system as the number of free parameters is far larger than 
the number of observational constraints (especially given that 
the observed 21cm profile can be reproduced with a single disk 
galaxy model). High resolution radio interferometric observa- 
tions at 805 MHz will help refine the kinematic models and 
shed more light on the nature of the intervening lens system. 
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